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5 . Abstract 



The phase diagram of a two-fluid bosonic system is investigated. The proton- 
neutron interacting boson model (IBM-2) possesses a rich phase structure involving 
I three control parameters and multiple order parameters. The surfaces of quantum 

phase transition between spherical, axially-symmetric deformed, and SU*^(3) tri- 
axial phases are determined, and the evolution of classical equilibrium properties 
\ across these transitions is investigated. Spectroscopic observables are considered in 

■ relation to the phase diagram. 
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C ' 1 Introduction 



X 

^ I The phase structure of quantum many-body systems has in recent years been 

a subject of great experimental and theoretical interest. Models based upon 
algebraic Hamiltonians are well-suited to the study of phase transitions. They 
possess a well-defined classical limit [1], allowing classical order parameters 
to be determined. And for certain specific forms of their Hamiltonians, alge- 
braic models exhibit dynamical symmetries, which correspond to qualitatively 
distinct ground-state equilibrium configurations. These constitute the phases 
of the system [2]. Algebraic models have found extensive application to the 
spectroscopy of many-body systems, including nuclei [3] and molecules [4]. 

In the present work, the phase structure of a system comprised of two inter- 
acting fiuids is investigated. The phase structure of one-fiuid algebraic models, 
especially the interacting boson model (IBM) [3] for nuclei, has been studied in 
detail [2, 5]. However, algebraic models may also be used to describe multi-fluid 
systems, with multiple interacting constituent species. While one-fluid systems 
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are described by a single elementary Lie algebra, usually U(n), multi- fluid sys- 
tems are described by a coupling of such Lie algebras, Ui(n)cg)U2(^)(S'- • ■ [3,4]. 
A richer phase structure arises for multi- fluid models, not just from the greater 
number of control and order parameters afforded by the more complicated 
model, but more fundamentally from the coupling of multiple subsystems, 
each of which can exist in a different phase or can drive phase transitions of 
the composite system. 

Here we consider the proton- neutron interacting boson model (IBM-2) [3, 6-8], 
in which proton pairs and neutron pairs are treated as distinct constituents. 
The one-fluid IBM, with U(6) algebraic structure, has three dynamical sym- 
metries, separated by flrst and second order phase transitions [2,5]. The IBM- 
2, with U(6) ® U(6) algebraic structure, supports four dynamical symme- 
tries [9,10]. Its phase diagram is therefore more involved and is found to 
possess qualitatively new features. Due to the complexity of the problem, 
a combination of analytic and numerical methods have been applied in the 
present work. Preliminary results were reported in Refs. [11,12]. Numerical 
studies of the IBM-2 phase structure have also been carried out by Arias, 
Dukelsky, and Garcia-Ramos [13,14]. 

The IBM-2 and its classical limit are summarized in Sec. 2. The phase diagram 
of the IBM-2 is flrst investigated for an essential Hamiltonian with few pa- 
rameters, for which the most complete analytic results are obtained (Sec. 3). 
The treatment is then extended to a more general Hamiltonian, incorporat- 
ing realistic quadrupole and Majorana interactions (Sec. 4). Connection of 
the IBM-2 phase diagram with experimental data requires knowledge of the 
spectroscopic predictions across the phase transitions. These are addressed in 
Sec. 5. 

The results of Sec. 5 are of relevance in the search for triaxial shapes in nu- 
clei. Speciflc signatures of two-fluid triaxial deformation, and of the phase 
transition between axially symmetric and triaxial structure, are presented. 
Those signatures involving low-lying states are applicable both to current ex- 
periments and to experiments planned for next generation radioactive beam 
facilities. Those involving the high-lying magnetic dipole mode might be most 
directly investigated through resonance fluorescence experiments at high in- 
tensity gamma-ray facilities. 
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2 IBM-2 definition and classical limit 

2. 1 Hamiltonian 

Let us first summarize the IBM-2 Hamiltonian and the dynamical symmetries 
it supports. Operators in the IBM-2 are constructed from the generators of 
the group U7r(6) (8) Uj/(6), realized in terms of the boson creation operators Sp o 
and d'jj^^ (where p represents tt or u, and p=—2, .... 2) and their associated 
annihilation operators, acting on a basis of good boson numbers N^^ and N^,. 
The physically dominant interactions are contained in a Hamiltonian 

H = e^fld^ + e^Udu + «,r7r<5f • + «7r.Qf " Q^" + ^^vvQl" ' Q^" + AM, (2.1) 

where nap = 4 ■ dp, = {sj, x + d^ x 5^,) (2) + ^pid], x dp)^^\ M = 

conventional spherical tensor coupling notation {e.g., Ref. [15]) has been used. 
This Hamiltonian contains one-fluid contributions, arising from like-nucleon 
pairing (ridp) and quadrupole {Qp ■ Qp) interactions, as well as two-fluid cou- 
pling terms, arising from proton-neutron quadrupole (Qtt • Qu) and Majorana 
(M) interactions. The physically relevant ranges of the Hamiltonian parame- 
ters are £p > 0, Kpp> < 0, -a/7/2 <Xp< ^/7/2, and A > [3]. 

A dynamical symmetry occurs when, for certain values of the parameters, 
the Hamiltonian is constructed from the Casimir operators of a chain of sub- 
algebras of U^(6) ® Ui/(6). Three of the IBM-2 dynamical symmetries oc- 
cur for Xtt — Xv and have direct analogues in the one- fluid IBM [9]. When 
/^TTTT = ^^■Ku = i^vv — 0, the U7ri/(5) dynamical symmetry is realized, with subal- 
gebra chain 

U^(6) ® U,(6) D U^(5) (8) U,(5) D U^.(5) D S0^,(3). (2.2) 

The geometric interpretation is that the proton and neutron fluids undergo 
oscillations about a spherical equilibrium conflguration. Parameter values £^ = 
Sn — Q with Xtt — Xv — ^ produce the S0^i,(6) dynamical symmetry 

U^(6) ® U,(6) D S0^(6) ® S0,(6) D S0^,(6) D S0^,(3), (2.3) 

yielding deformed, 7-unstable structure. And — e„ — Q with Xtt — Xv — 
— -\/7/2 gives the SU7r;/(3) dynamical symmetry 

U^(6) ® U,(6) D SU,(3) ® SU,(3) D SU,,(3) D S0,,(3), (2.4) 

for which prolate axially symmetric structure is obtained. The complemen- 
tary case Xtt = Xi/ = +\/7/2, giving oblate axially symmetric structure, is 
distinguished by the notation SU7ri/(3). 
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Fig. 1. Geometrical interpretation of the equilibrium structure for the SU*^(3) dy- 
namical symmetry. A prolate deformed proton fluid (light gray) and oblate de- 
formed neutron fluid (dark gray) are coupled with orthogonal symmetry axes. Figure 
adapted from Ref. [22]. 

However, a symmetry special to the IBM-2, denoted SU*j,(3), is obtained for 
= = with Xtx = — \/7/2 and Xi^ = +\/7/2 [10, 16-19]. In this case the 
Hamiltonian is constructed from Casimir operators of the subalgebra chain 

U^(6) U,(6) D SU^(3) SUA^ D SU;,(3) D S0^,(3). (2.5) 

The equilibrium configuration consists of a proton fluid with axially symmetric 
prolate deformation coupled to a neutron fluid with axially symmetric oblate 
deformation, with their symmetry axes orthogonal to each other [10,20,21]. 
This yields an overall composite nuclear shape with triaxial deformation, as 
shown in Fig. 1. To avoid ambiguity, we shall adopt here the notation 811*^,(3) 
for the complementary case x-k — +'\/7/2 and Xu — /2, in which the proton 
and neutron deformations are interchanged. 

2.2 Coherent state energy surface 

The phase structure of an algebraic model is determined from the classical 
hmit of the model, which is obtained through the coherent state formal- 
ism [1,23,24]. A coherent state of good boson number is constructed, and 
the parameters of this state are directly related to the classical coordinates of 
the system. This coherent state is used as a variational trial state: the global 
minimum of the energy surface S = {H) is used to determine the ground 
state energy and equilibrium classical coordinates of the system. In the limit 
of infinite boson number, the coherent state variational energy converges to 
the true ground state energy [1,23]. Detailed coherent state analyses of the 
one-fluid IBM may be found in Refs. [2,5,25-30]. 
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Once the equilibrium properties of the model have been established, ground 
state phase transitions can be identified and categorized according to the 
Ehrenfest classification [31]: a phase transition is first order if the first deriva- 
tive of the system's energy is discontinuous with respect to the control pa- 
rameter being varied, second order if the second derivative is discontinuous, 
etc. Where the system's energy is obtained, as in the present classical anal- 
ysis, as the global minimum of an energy function S, a first order transition 
is usually associated with a discontinuous jump in the equilibrium coordi- 
nates ("order parameters") between distinct competing minima. Second or 
higher order transitions are associated instead with a continuous evolution of 
the equilibrium coordinates, as when an initially solitary global minimum be- 
comes unstable (possessing a vanishing second derivative with respect to some 
coordinate) and evolves into two or more minima. It should be noted that, 
whenever the order of a phase transition is obtained by numerical analysis, 
application of the Ehrenfest criterion is limited by the ability to numerically 
resolve sufficiently small discontinuities. This is especially a consideration for 
points of first-order transition very close to a point of second-order transi- 
tion. Moreover, problems with the classification scheme, not addressed here, 
may arise at the boundaries of the parameter space or when the Hamiltonian 
possesses additional symmetries. 

For a two-fiuid model, the boson number for each constituent fiuid is conserved 
separately. The coherent state is constructed not only with good total boson 
number but with good boson number for each fluid. For the IBM-2, the number 
operators are Np = o'5p,o + S^=-2 ^\,ij.^p,ij.- The IBM-2 coherent state [26, 32- 
34] is deflned in terms of proton and neutron condensate bosons 

where |af | = (EJ=-2 

lAT^, 4^); iV., . {Blf^ K)"^ |0). (2.7) 

The a^pji^ are interpreted geometrically as quadrupole shape variables [35] for 
the proton and neutron fluids. They are related to the four deformation pa- 
rameters (/Stt, Ttt, Pu, and j^) and to the six Euler angles {Oitt, 02^, ^stt, 0iv, 
62U, and 6^u) specifying the orientations of the proton and neutron intrinsic 
frames by 



+ -^/5,sin7,[D2*(^i„^2p,M +^J-2(^ip,^2p,Mj> (2-8) 



where -D^/m is the Wigner D function [36] . 
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Calculation of the expectation value of H with respect to the coherent 
state (2.7) yields an energy surface S = {N^, a^^); a'^^'^\H\N^, a^^); N^,, a^^'^). 
By rotational invariance, £ can only depend upon the relative Euler angles 
between the proton and neutron fluid intrinsic frames, not the and 9ii, 
separately. For evaluation of £, the Euler angles may thus simply be chosen 
to be 6i„ = and 9ii, — di. Investigations of the IBM-2 coherent state energy 
surface have been carried out in Refs. [20, 21, 37]. 

The expectation value with respect to the IBM-2 coherent state of an operator 
constructed from the boson operators of one fluid only {e.g., hdn or ■ Qtt) 
can be calculated as in the one-fluid IBM, by the methods of Refs. [26,28]. 
The expectation values of the one-fluid operators appearing in (2.1) are [28] 



5 + (1 + xM 



(2.9) 
(2.10) 



+ 



NpjN, - 1) 

(1 + Pir 



4/5p'-4W-Xp/5,^cos37, 
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The expectation value of a two- fluid operator constructed as the scalar product 
of one-fluid operators {e.g., Qtt • Qv) can be obtained using a factorization 
result [21, (CI)]. However, for more complicated operators, the general method 
presented in Appendix A provides a convenient means of calculation. The 
expectation value of Q.„ -Q^ is a. function of all seven possible coordinates (/^jr, 
Ttt, Pu, lu, A, 1^2, and T?3), 



[4^)* + 4') + x.(4^)*x4^))(^)] 

•[«t^^* + «l^^ + X.(4')*x4^))(^)], (2.11) 



with the a^^^ as deflned in Eqn. (2.8), which for vanishing Euler angles sim- 
phfles to 



7(l + /32)(l + /32) 

VuIx^Ptt cos(27^ + 7,,) + XuPu cos(7^ + 27^,)] 



14 cos(7^ - 7,.) + XttXuPttPu cos(277r - 27^) 

(2.12) 



The expectation value of the Majorana operator is a complicated function of 
the coordinates. For vanishing Euler angles, it reduces to 



(M) 



A^A^ 



{1 + Pl){l + P2) 



I3i + I3t- 213,13, cos(7, - 7,) + /?^/?^ sin^(7, - 7, 



(2.13) 
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A detailed study of the possible equilibrium Euler angle and 7^ values for the 
most general two-body lBM-2 Hamiltonian has been presented in Ref. [21]. So 
long as the multipole decomposition of the proton-neutron interaction contains 
no hexadecapole component [{dl x dj^Y'^^ ■ {dl x d^)'^^^], it is found that the equi- 
librium configuration must have Euler angles which are vanishing or multiples 
of 7r/2. Thus, the proton and neutron intrinsic frames are "aligned", to within 
a possible relabeling of axes. [For clarity, we note that in the SU*j,(3) equi- 
librium configuration (Fig. 1), even though the proton and neutron symmetry 
axes are orthogonal, the intrinsic frame coordinate axes are actually parallel, 
i.e., the Euler angles vanish. The proton (7^ = 0) symmetry axis is the z 
axis, while the neutron (7,^ = 7r/3) symmetry axis is the y axis.] If, instead, 
a hexadecapole contribution is present, "oblique" equilibrium configurations 
are possible. The operator " Qu contains no hexadecapole component, but 
the Majorana operator does (Sec. 4.3). The IBM-2 is thus characterized by ei- 
ther four order parameters (/3^, 7^, P^,, and 7,^) or seven order parameters (/3^, 
Ttt! I^u, 1u, '&2, and ■j^a), depending upon the form of the proton- neutron 
interaction. 



3 Phase structure for an essential Hamiltonian 

3.1 The F-spin invariant Hamiltonian 

A simple, schematic Hamiltonian which retains the essential dynamical sym- 
metry features of the IBM-2, obtained as a special case of Eqn. (2.1), is 

H = e{nd. + Udu) + k{Q^- + ) ■ (Qf + Q^)- (3-1) 

This Hamiltonian is often preferred for theoretical studies due to its invariance 
under F-spin rotations {e.g., Ref. [38]) for Xt^ = Xi^- To make the analysis of 
the IBM-2 phase structure as transparent as possible, let us first consider in 
detail this F-spin invariant Hamiltonian before proceeding in Section 4 to the 
more general form. It is also convenient for the moment to restrict attention 
to the case N^^ — N^,. 

To obtain the classical limit, it is necessary to take N ^ 00, where A^" = 
Nt^ + N^. In this limit, the lower-order terms with respect to A^ in Eqn. (2.10) 
become negligible, and it is seen that the quantities (n^p) are linear in A^, 
while the (Qp ■ Q^i) are all quadratic in A^. To prevent the {fidp) contributions 
from vanishing identically when the limit is taken, it is necessary to rescale 
the Hamiltonian coefficients by appropriate powers of A". It is thus convenient 
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su;;,(3) 



Fig. 2. Parameter space for the F-spin invariant Hamiltonian of Eqn. (3.2). 
For XV = (shaded plane) the analysis reduces to that for the one-fluid IBM 
(see text), yielding a point of second-order phase transition embedded in a 
curve (3.3) of first-order phase transition points. The axes cover ranges < ^' < 1, 
— \/7/2 < Xs ^ and — \/7/2 < xv < The xs a-^d Xv axes are scaled 

by so as to converge to a point at the V^ri/i^) limit, where the Hamiltonian is 
independent of xs and xv- Figure from Ref. [22]. 

to reparametrize the Hamiltonian (3.1) as 

H = ^{nan + nau) - ^{Q^^ + Ql") ■ (Qf + Q^l, (3.2) 

so that the energy function 8 is independent of A^. This definition also con- 
denses the full range of possible ratios e/k onto the finite interval < ^' < 1. 
An overall normalization parameter for H has been discarded as irrelevant to 
the extremum structure of the energy surface. 

There are three control parameters — x-ki ^^'^ Xv — for the Hamilto- 
nian (3.2). It is convenient to alternatively introduce "scalar" and "vector" 
parameters xs=\{x-k + Xv) ^'^'^ Xv=\{x-k ~ Xv) - The parameter space is out- 
lined in Fig. 2. 

For Xv — 0, indicated by the shaded triangle in Fig. 2, the analysis essentially 
reduces to that for the one-fiuid IBM. The equilibrium configuration occurs 
for jSj^ = f3u{= 13) and = 7jy = and is identical to that obtained for the IBM 
Hamiltonian iJiBM = [(1 — ^')/N]'f^d — {^'/N'^)Q^ ■ Q^, the phase structure 
of which is well known [2,5], with x = Xs- ^oi < 1/5, the energy surface 
has a minimum at (3 — 0. At ^' = 1/5, this minimum becomes unstable with 
respect to /3, i.e., d'^E/dp'^ = 0. If x = 0, this instabifity leads to a second- 
order transition between undeformcd (/? = 0) and deformed {j3 ^ 0) structure 
at ^' = 1/5. However, for any other value of X) the minimum at /9 = is 
preempted as global minimum by a distinct minimum with nonzero /? before 
^' — 1/5. This leads to a first-order phase transition, at a parameter value 



= i'cix) given by 

e(x) = (3.3) 

easily derived from Ref. [29, (2.9b)]. Thus, the point of second-order phase 
transition at ^' = 1/5 and x = hes on a trajectory (3.3) of first-order transi- 
tion points (Fig. 2). 

It simphfies the analysis of the phase diagram for the remainder of the param- 
eter space to note the presence of reflection symmetries with respect to the Xp 
parameters. Each term contributing to £ is seen by inspection of Eqns. (2.9), 
(2.10), and (2.12) to be invariant under a simultaneous transformation of co- 
ordinates and Hamiltonian parameters 

TT TT 

Thus, the equilibrium deformation at a point {xn, Xu) in parameter space is 
simply related to that at the point (— Xtt, —Xiy) by Eqn. (3.4). The phase dia- 
gram is symmetric under simultaneous negation of Xw and Xu o^) equivalently, 
simultaneous negation of xs and xv- Moreover, for the F-spin invariant Hamil- 
tonian (3.2) with A^TT — N^, the numerical coefficients on corresponding proton 
and neutron terms in £ are equal. Thus, £ is symmetric under interchange of 
all proton and neutron variables, and the phase diagram is symmetric under 
interchange of Xtt and Xiy or, equivalently, negation of xv- 



3.2 The SU^^(3) - S0^^(6) - SU;^(3) plane 



Let us begin the investigation of the phase structure of the Hamiltonian (3.2) 
with an analytic study for ^' = 1, corresponding to the rightmost plane of 
the parameter space diagram in Fig. 2. This plane encompasses the S07ri/(6), 
SU7ri/(3), and SU*^(3) dynamical symmetries. 

Surrounding the SU^,^(3) dynamical symmetry is a region of parameter space 
in which the equilibrium deformations are axially symmetric (7,^ = 7^ = 0), and 
a similar region surrounds the SU,r!y(3) dynamical symmetry (7,1 = 7;/ = 7r/3). 
In either of these cases (77r = 71/ = or 7r/3), the energy surface can be written 
in an especially simple form in terms of the function F defined in Eqn. (B.l), 



£:(/3,,/3,) = -4 



1 



F - 



1 



n 2 



(3.5) 



where a = cos37p = ±1. The energy function in (3.5) is minimized when 
the magnitude of the quantity within brackets is maximized. In the vicinity 
of SU^,^(3), including in the entire quadrant for which Xtv and Xi/ are both 
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negative, it follows from the results of Appendix B that the minimum in 
S{(3t^,Pi,) occurs at coordinate values 



Xp 



+ 1 



1/2 



(3.6) 



with a — +1 (i.e., = 7i/ = 0). Similary, in the vicinity of S\Jt^„{3), including 
in the entire quadrant for which Xn and Xv ^-re both positive, the minimum in 
^{Pir,Pu) occurs at coordinate values 



Pp 



Xp 



14 



Xp 



14 



+ 1 



1/2 



(3.7) 



with (7 = —1 {i.e., 7,^ = 71, = 7r/3). In either case, the value of the energy 
function at the global minimum is 



with and Pi, given by Eqn. (3.6) or (3.7) as appropriate. 



(3.^ 



As x-K a.nd Xu a-^s varied away from their SU7ri/(3) values, eventually axial 
equilibrium deformation gives way to triaxial deformation, with 7^ and/or 7^ 
nonzero. This transition occurs on the locus of points at which the minimum 
given by (3.6) first becomes unstable with respect to 7 deformation. Since S 
depends upon both 7v and 7j,, instablility occurs when the directional second 
derivative of S first vanishes along some "direction" in {'jTr,lu) coordinate 
space, which may generally happen before either d'^S/dj^ or d'^S/d^l vanishes 
individually. The equation describing the boundary curve in Xn and Xu is most 
compactly expressed in terms of the corresponding equilibrium values and 
Pi, from Eqn. (3.6) or (3.7), as 



33 fl3 



(3.9) 



This curve is shown in Fig. 3. The direction in (77r,7!y) space in which the 
instability occurs, at any given point on this boundary, is given by 



where 



7. ^ l + Pl B ^ \-rPl 2C 
7, \ + Pl2A I + PIB' 



A ^ gplipl - 1) + PnPA'^Pl - iM + 1) 

B^-Af3MP'-^M-^) 

c ^ 9pl{pl - 1) + PM^pl - iM + !)• 



(3.10) 



(3.11) 



Along the SU^i,(3)-SU*^(3) line in particular, for which x-k = — \/7/2, the 
transition occurs at Xi' ~ 0.4035, obtained from Eqn. (3.9) as the root of a 
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Fig. 3. Phase diagram for the SU7ri/(3)-S07ri/(6)-SU*y(3) plane = 1) in the pa- 
rameter space of the Hamiltonian of (3.2), for A?^ = showing the curve (3.9) of 
second-order phase transition. The diagram is rotated to ahow more direct compar- 
ison with Figs. 2 and 4. Figure adapted from Ref. [12]. 

quartic equation. In this special case, the global minimum becomes soft purely 
with respect to at fixed = 0. 

The analytic results of this section indicate a continuous transition between 
axially symmetric and triaxial structure on the curve (3.9). In Sec. 3.4 it is 
verified numerically that this transition is second, rather than higher, order, 
and that no prior phase transition to a distinct minimum occurs before the 
curve (3.9) is reached. 

3.3 The transition between spherical and deformed equilibrium 

The transition between spherical and deformed equilibrium for the IBM- 
2 for arbitrary xv is closely related to the transition which occurs in 
the special plane Xv = (shaded in Fig. 2). The energy function 
^ (C', X5, Xv\ Pit, Ttt, lu) obeys the identity 

xs, XV, 0, /9, 0) = S{e, xs, 0; /3, 0, /3, 0), (3.12) 

i.e., for a deformation which is axial and has = /3^, the value of £ is 
independent of Xv- Recall that the equilibrium deformation for xv = is 
always axial with /5,r = (Sec. 3.1). It is straight-forward to demonstrate from 
these observations (see below) that a nonzero equilibrium deformation for the 
parameter point (^', xs) in the Xv — plane implies a nonzero equilibrium 
deformation for all the out-of-plane points {^',XS:Xv)- Thus, the ^' value 
C'cixs, Xv) for which the spherical-deformed transition occurs at a given xs 
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and xv obeys, by Eqn. (3.3), 



C(xs,Xv) < 



5 + fxl 



(3.13) 



The detailed argument proceeds as follows. Let the equihbium energy at a gen- 
eral point in parameter space be denoted by So{^' , xs , Xv) , and let Po{^',Xs) 
be the equilibrium value of (3-,^ and [3^ for xv = as determined, e.g., from 
Ref. [29, (2.12)]. Then SoiC, Xs, 0) =£[e, Xs, 0; PoiC, Xs), 0, /3o(e', Xs), 0]. The 
global minimum energy for arbritrary xv then is subject, by Eqn. (3.12), to 
the upper bound 

^o(e', XS: Xv) < XS, XV, m', xs), 0, m', xs), 0] . 

^S[e,xs,0;f3o{e,xs),0,f3o{e,xs),0] = £o{e,xs,0). ^ 

Note that S evaluated at zero deformation is zero for the Hamiltonian con- 
sidered here (2.1), so < if and only if the configuration is deformed. If 
/3oi^',Xs) 0, then £oi^',Xs,Xv) < for all xv by Eqn. (3.14), and the 
equilibrium deformation for parameter point {i',Xs,Xv) is also nonzero. 

To further investigate the nature of the spherical-deformed transition, let us 
consider the stability of the minimum of S at i3t^ = f3,^ = 0. Instability occurs 
when the second derivative of 8 vanishes along some "direction", which we 
denote (ut^jU^), in (Pt^jP^) coordinate space, for some fixed values of and 
7,^. The second derivative of £ along the ray parametrized as = u^^P and 
= Ui,f3 is 

= (1 - 30{ul + ul) - 4C'unU,cos{j^ - 7,). (3.15) 

/3=0 

This is independent of xs and xv and depends upon the 7p coordinates only 
through their difference — ■y^. The smallest ^' value at which the second 
derivative (3.15) vanishes is ^' = 1/5, with — Ui,, indicating instability of the 
minimum against deformations with (3-^ — (3,^ and 7^ = 7i/(= 7)- As the second 
derivative is independent of 7, the system is unstable against deformations 
of all possible 7 values, representing prolate, oblate, and intermediate triaxial 
deformations, simultaneously. 

For X5 7^0, it is apparent from Eqn. (3.13) that the spherical-deformed phase 
transition occurs before the spherical minimum becomes unstable (^' = 1/5), 
and so the transition is a first-order transition to a distinct minimum. Thus, 
the curve of first-order phase transition occuring in the plane Xy = "prop- 
agates" out of this plane to form a surface of first-order transition. In the 
approximation that /^^r ~ f3u for the deformed minimum, the surface would be 
a vertical extension of the one-fluid IBM transition trajectory in Fig. 2. The 
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deviation of the surface from this hmiting location is studied numerically in 
Sec. 3.4. 

For xs — 0, however, it is possible for the minimum at zero deformation to 
remain the global minimum until ^' = 1/5, leading to a second-order phase 
transition. Let us consider the hkely properties of any alternative lower, de- 
formed global minimum (these arc confirmed numerically in Sec. 3.4). By 
symmetry in the proton and neutron parameter magnitudes {x^k = —Xv)-, it is 
reasonable to expect this minimum to have (3^^ — (3u- A minimum with axial 
deformation has been excluded by Eqn. (3.3), but a minimum with triaxial 
deformation must be considered. Since Xs = ^ represents the "reflection" plane 
in parameter space between prolate and oblate structure (3.4), it is natural to 
expect the minimum to have (77r + 7;/)/2 = ti/6. Inspection of the restricted 
energy function S{^',0,xv', f^j'K/fi + 7i/,/3, 7r/6 — jv) reveals that a triaxial 
deformed minimum docs preempt the minimum at zero deformation as global 
minimum [at ^^(0, Xv) = V (3 + 7x^^/2-1-2x^/7)] but that this only occurs for 
Ixvl ^ ^7/2, which lies outside the parameter range |xy| ^ '\/7/2 of interest 
in the present study. 

In conclusion, a line of second-order phase transition between spherical and 
deformed equilibrium occurs at ^' = 1/5 and xs = 0, subsuming the one- 
fluid IBM second-order transition point. This line is embedded in a surface of 
first-order transition points. 

3.4 Numerical investigation of the full phase diagram 

The remainder of the phase diagram is obtained by numerical minimization 
of the energy surface with respect to [3.,^, 7,^5 Pu, and 'y^. This minimization 
provides the equilibrium coordinate values at any point in the parameter space 
of Fig. 2. 

For robust identification of the global minimum, S is first evaluated at each 
point on a fine mesh in the deformation coordinates {Aj3p = 0.02 and A'jp = 1° 
provides sufficient resolution for most of the calculations shown here). All 
points which are discrete local minima of £ relative to the neighboring mesh 
points are identified. The coordinate values for these minima are then refined 
iteratively by a conjugate gradient method. The global minimum is identified 
from among the refined values. 

The phase diagram for the F-spin invariant Hamiltonian with N^^ — N^, ob- 
tained numerically in this fashion, is shown in Fig. 4. [Only one quadrant of 
parameter space is included in this plot, since the others may be obtained 
by refiection (Sec. 3.1).] The line of second-order transition between spherical 
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Fig. 4. Phase diagram of the proton-neutron interacting boson model (IBM-2) for 
the F-spin invariant Hamiltonian (3.2) with iV^ = N^,, as obtained by numerical 
minimization of £. The surfaces of first-order and second-order transition between 
regions of undcformcd, axially symmetric deformed, and triaxially deformed equilib- 
ria are shown. Only one "quadrant" of parameter space (0 < ^' < 1, —V7/2 < xs ^ 0; 
and < xv < V^/^) is included in this plot, since the others may be obtained by 
reflection. The xs and xv axes are scaled by so as to converge to a point at the 
U^,y(5) hmit. Figure from Ref. [12]. 

and deformed equilibrium at ^' = 1/5 and X5 = is the locus of simultaneous 
contact of four regions of the phase diagram, those of spherical, prolate axi- 
ally symmetric, oblate axially symmetric (not shown, in adjacent quadrant), 
and triaxial equilibria. Numerical investigation of the behavior of £ at points 
on the axial-triaxial transition surface allows the Ehrenfest criterion to be 
applied, and it appears that the transition is everywhere second order. 

We now consider in detail the evolution of the equilibrium energy and coor- 
dinates at selected locations throughout the phase diagram. Within the plane 
^' = 1, the boundary curve separating axially symmetric and triaxial equilib- 
ria was established in Sec. 3.2. The equilibrium properties along the 811^,^(3)- 
SU*j,(3) transition line are shown in Fig. 5. Several characteristics may be 
noted. The second derivative of £ with respect to the control parameter Xu 
is discontinuous, as in an Ehrenfest second-order phase transition. The first 
derivative of one order parameter, 71,, is infinite at the critical point, with an 
approximately square-root dependence upon the control parameter after this 
point, as in a Landau second-order phase transition [39]. The first derivative 
of at this point is, in contrast, discontinous but finite. In the special case 
of the SU7ri/(3)-SU*^(3) fine, as already discussed, the two other order param- 
eters Pt, and completely decouple from the transition process, remaining 
constant throughout. 

In Fig. 6, the equilibrium properties on lines of constant xv in the plane 
— i, which would appear as horizontal lines across the diagram of Fig. 3, are 
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Xv 

Fig. 5. Evolution of the equilibrium values of (a) £, (b) {3^-, and (c) '^^ along the 
line in parameter space between the SU7ri.(3) and SU*,^(3) dynamical symmetries 
(Xtt = — — \/7/2 <Xy^ obtained numerically. The dashed line indi- 

cates the second-order transition point, Xv ~ 0.4035. 

shown. These hnes cross both the prolate-triaxial and triaxial-oblate second- 
order transition curves. The equihbrium coordinate values for positive xs ^-re 
related to those for negative xs by the transformations of Sec. 3.1. Observe 
that the order parameters (/Stt, Ttt, f^v, and 7,^) are coupled in their behavior 
at the transition points, all simultaneously exhibiting discontinuities in their 
derivatives at each transition point. In the limiting case Xy = [Fig. 6(right)], 
the two points of second order transition at positive and negative xs converge 
to a single point at xs = 0. The discontinuities in the second derivative of £ 
with respect to xs at the two transition points [Fig. 6(left, middle)] combine 
to form a "cusp" in S. The discontinuity oi d£/dx at x = 'w^as noted in 
the context of the one-fluid IBM by Jolie et al. [40] . Such a discontinuity in 
the slope of £ would, according to the Ehrenfest criterion, indicate a first- 
order phase transition. However, the present case provides an example of the 
great variety of phenomena possible in multi-parameter problems, as is does 
not naturally fit the one- dimensional Laundau model for a first-order phase 
transition. The conventional Landau first-order phase transition [39] involves 
"competition" between two distinct local minima simultantously present in 
the energy surface, one overtaking the other as global minimum; whereas, 
in the present example, a single minimum with 7,^ = 7jy = is present for 
X5 < 0, a single minimum with 77r = 7^ = vr/3 is present for xs > 0, and a 
continuous trajectory of minima of arbitary 7(= 7^ = ji,), arising from the 
S0,ri/(5) symmetry of the Hamiltonian, are present for xs — 0. This example 
suggests the need for a comprehensive extension of the classification schemes 
for single-parameter phase transitions to cover phase transitions with multiple 
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Fig. 6. Evolution of the equilibrium properties along rays of constant xv in the =1 
plane of parameter space, passing from the region of axially symmetric prolate equi- 
librium to that of axially symmetric oblate equilibrium. Results are shown for the 
F-spin invariant Hamiltonian (3.2) with N-,j = Ny, along the rays of (left) = 0.8, 
(center) xv = 0.4, and (right) xv = 0. Proton fluid variables (/^tt and 7,^) are rep- 
resented by solid curves, while neutron fluid variables {j3y and 7,^) are represented 
by dashed curves. The two points of second order phase transition (left, middle), 
converging towards the single point X5 = in the limit xv = (right), are marked 
with arrows. 

order parameters (see also Ref. [41, Chapter 17]). 



Finally, we consider the transition between spherical and deformed equilibrium 
configurations. In Sec. 3.3, an upper bound, i'^ixsiXv) ^^'cixsi^) was placed 
upon the location of the spherical-deformed transition surface [see Eqn. (3.13)], 
equality holding in the case (3.,^ = (3^. The actual shape of this surface, visible 
on a coarse scale in Fig. 4, is shown in detail in Fig. 7. Several "vertical" slices 
through the surface, at constant xs, are plotted. The hue of second-order 
phase transition at = 1/5 appears rightmost in the figure. The curves vary 
little in ^' (< 0.2%) below their upper bound. 

The evolution of equilibrium properties across the spherical-deformed transi- 
tion is shown in Fig. 8. The second-order transition along the U7ri^(5)-SO^,^(6) 
line and the first-order transition along the U7ri/(5)-SU7ri^(3) line, quantitatively 
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Fig. 7. Slices through the boundary surface between the regions of spherical and 
deformed equilibrium at fixed xs- The ^' value at xv = is the value ^'dxs) from 
the one-fluid IBM (3.3), and the recession of each curve towards smaller ^' above 
Xy = indicates the extent to which ^^(xs, Xv) recedes below the upper limit (3.13) 
obtained using (3-,^ = The curves shown are for X5 = to 1 in steps of 0.1 (right 
to left). 

identical to the corresponding transitions in the one-fluid IBM, are shown for 
reference [Fig. 8(left, middle)]. The transition along the U7rj,(5)-SU*j,(3) line 
[Fig. 8 (right)] is strongly constrained by the proton-neutron interchange sym- 
metry about X5 = 0, with /J^r = and with 7^ and 71/ symmetric to each other 
about 7r/6. 



4 Phase structure for a general Hamiltonian 



^.1 Energy surface parameters 



In the preceding section, the general features of the IBM-2 phase diagram were 
estabhshed using a schematic Hamiltonian, and only the case — N^, was 
considered. We now address the more general Hamiltonian of Eqn. (2.1) and 
also consider arbitrary values of the ratio N-,^/Ny. These two generalizations 
are closely related, as the coefficients in the energy surface £ depend upon 
the Hamiltonian coefficients and the boson numbers in combination. To make 
explicit the dependence of each term contributing to E upon parameters, boson 
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Fig. 8. Evolution of tlic equilibrium properties in the vicinity of the transition be- 
tween spherical and deformed equilibria. Results are shown for the F-spin invariant 
Hamiltonian (3.2) with A^^ = A^j,, along (left) the U-n-i/{5) — S0.^i,{6) transition line 
ixs = 0, XV = 0), (middle) the U7rzy(5) — SU7rjy(3) transition line {xs = —V7/2, 
XV = 0), and (right) the \Jtti/{5) — SU*^(3) transition line {xs = 0, xv = —V7/2). 
Proton fluid variables are represented by solid curves, while neutron fluid variables 
are represented by dashed curves. Equilibrium values are undefined, and there- 
fore not shown, when Pp = (h,i) and along the entire U7ri/(5) — S07ri/(6) transition 
line (g). Points of phase transition are marked with arrows. 

numbers, and coordinates, we express £ as 

+ K^uN^N^ fnuiX^, Xu, Pw, 7^, Pu, 7^, ^1, ^2, ^3) + KuuN^ fuuiXu] Pu, 7u) 

+ XN^K fuiP.: 7^, 7., ^1:^2: ^s)- (4.1) 

The functions / {e.g., = (ndir) /N.,^) may be read off directly from 
Eqns. (2.9)-(2.13). The limit of both large N,^ and large iVj, has been taken, so 
that in {Qp ■ Qp) the term hnear in Np is suppressed. The energy surface (4.1) 
is entirely determined by the boson-number-weighted coefficients e'^ = SpNp, 
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K'pp, = Kpp'NpNpi, and A' = XNt^Ni,, together with the Xp, in terms of which 

S = e'J^{(3^) + e'J'APu) + nl^UAXn, 13^, In) 

Xu] Pn, Piy, 7v, ^^1, ^^2, '&3) + f^'yufvv{Xv\ 13 

+ A7m(/3.,7.,/5.,7.,^i,^2,^?3). (4.2) 
4..2 Quadrupole interaction coupling coefficients 

The dominant role of the proton-neutron quadrupole interaction in producing 
collective nuclear deformation has been well established [42]. Consequently, the 
like-nucleon quadrupole interactions are often neglected in the IBM-2 Hamil- 
tonian (see Refs. [7, 38, 43]). Microsopic estimates suggest that the shell model 
proton-proton and neutron-neutron quadrupole interactions each have ~1/10 
to 1/5 the strength of the proton-neutron quadrupole interaction [44]. How- 
ever, within the IBM-2, significant further strength is added to the effective 
like-nucleon interactions by renormalization effects arising from elimination of 
gi-wave nucleon pairs from the model space. This may yield like-nucleon cou- 
pling strengths comparable to the proton-neutron coupling strength [45] . The 
actual coupling strengths are a subject for further phenomenological study. 
Examples spanning a considerable range of values for Kpp/ are included in 
the following analysis. 

For investigation of the phase diagram for the Hamiltonian (2.1), it is conve- 
nient to again define a transition parameter ^' controling the relative weights 
of the hd operator and quadrupole operator in the Hamiltonian. This provides 
a coordinate system for the parameter space like that in Fig. 2. The contribu- 
tions of the different terms within each of these operators can then be specified 
by parameters and k'^p, defined such that 

yielding a A^p-independent expression 

E^{\- O {e'J^ + e'M - e + + KM (4.4) 

for the energy surface. To unambiguously define ^' and avoid redundancy 
in the parameters, we adopt the normalization convention + e'^, = 1 and 
Kn + Ku + Ku = ^- This choice gives e'^ = e'^j (e^ + e'^) and k'^^, = K,'pp, / {k,'^^ + 
k^^+k'^^) and is consistent with the definition of ^' in Eqn. (3.2). For the F-spin 
invariant Hamiltonian with arbitrary Nt^/N^, the parameters are e'^ = N^^/N , 
e', = NJN, k'^^ = Nl/N\ C = 'iN,NjN\ and k',, = Nl/N\ 



Ql'^]. (4.3) 
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Analytic results for the phase structure of the Hamiltonian (2.1) with arbi- 
trary quadrupolc coupling coefficients can be obtained in the case = 1 by 
a straight- forward extension of the analysis described in Sec. 3.2. Surround- 
ing the SU7„/(3) and SU7ri/(3) dynamical symmetries are regions of parameter 
space in which the equilibrium deformation is axially symmetric {'j.,^ — — 
or 7r/3). The energy surface in this case can be expressed in terms of the 
function F (of Appendix B) as 



, (4.5) 

with a = cos37p = ±1. This expression is simply a quadratic form in the 
quantities F, and the extremization problem reduces to a constrained mini- 
mization of the quadratic form with respect to these quantities. Since the Pp 
are limited to positive values, the F are restricted to the rectangular region < 
F{-Xn/VU;P.) < x+i-xjVTi) and < F{-Xu/VU; P,) < x+{-xjVU) 
for the case a = +1 [SU^j.(3)-like] or to x_{-x^/VU) < F{ -Xn/\^ ; -f^n) < 
and x_{~xjVl4:) < F(-x^/v^; -/5^) < for a = -1 [SU^^(3)-like]. Pro- 
vided Ktttt, Uttu, and K,i,u are all negative, the global minimum is again given by 
Eqn. (3.6) for a — +1 or Eqn. (3.7) for a — —1. The energy at the minimum 
is, in terms of the equilibrium Pp values, 

S = k'^X + k'^M + <Jl (4.6) 



The boundary curve separating axial and triaxial deformations is found as 
in Sec. 3.2 from analysis of the second derivatives of £. The axial minimum 
becomes unstable with respect to 7 deformation along a curve in parameter 
space, again most conveniently expressed in terms of the equilibrium Pp values, 
given by 

is^PM - 1) + M^Pl - + 1) 

2/3,(/32 - 2)2 

18^/5.(/52-l)+/3.(2/32-l)(/32 + i) 

X — — . (4.7) 

2/5.(/32-2)2 ^ ' 

The instability occurs against deformations with a ratio jn/ju given by 
Eqn. (3.10) with the values 

A ^ kjm - 1) + PM^Pi - ml + 1) 

B ^ -2K:j^(3,{(37r' - 2){Pl - 2) (4.8) 
C ^ 9kIMI - 1) + /3./?.(2/32 - l)(/?2 + 1). 

Along the SU7ri/(3)-SU*y(3) line in parameter space, the location of the phase 
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Fig. 9. The curve of second-order phase transition between axial and triaxial de- 
formation in the plane ^' = 1, shown for the i^-spin invariant Hamiltonian with 
different boson number ratios Nt^/Ni,: (a) N.,^/Ny = 2 and (b) N-j^/Ny = 4. The 
boundary curve for NT^/Ny = 1 (Fig. 3) is shown in each panel (dashed curve) for 
comparison. 

transition is obtained by solution of the quartic equation 

IS'^WI - 1) + - ml + 1) = 0, (4.9) 

and a similar relation with interchanged proton and neutron labels holds for 
the SU7riy(3)-SU*^(3) line. Observe that the value at which the transition 
occurs depends only upon the ratio k^^^/k^^. 

The dependence of the second-order transition curve upon Nt^/N^ is illustrated 
in Fig. 9. Since ^ k'^^^^ the phase diagram is not symmetric under reflec- 
tion with respect to xs or xv but remains symmetric under inversion (3.4). It 
is seen that increasing N^/N^ from unity moves the SU^,^(3)-SU*j,(3) transi- 
tion farther away in parameter space from 811^^(3) while moving the SU^iy(3)- 
SU*j,(3) transition closer. The Nt^/N^ dcpcnccncc of the boundary curve raises 
the possibility of considering phase transitions at a fixed point in the Hamilto- 
nian parameter space with N^^/Ny as the control parameter, i.e., phase tran- 
sitions as a function of particle number. 

For the F-spin invariant Hamiltonian considered in Sec. 3, the transition be- 
tween prolate axially symmetric equilibrium [SU7r;y(3)-likc] and oblate axially 
symmetric equilibrium [SU^,^(3)-like] always proceeds through an intermedi- 
ate stage of triaxial equilibrium, except at the point in parameter space cor- 
responding to S07riy(6). However, if k'^^ and sufficiently small relative 
to K^^, it is possible for a direct transition to occur between the prolate and 
oblate equilibria of Eqns. (3.6) and (3.7) without cither minimum becoming 
unstable with respect to triaxiality and thus for the regions of prolate and 
oblate equilibrium in the phase diagram to share a common boundary. The 
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boundary curve is given by 



/2 2/'2ii\ /2 2/2ii\i / / /2ii\ 

- <,<,a,a,(a2 + 1) - ^Hal - al) = 0, (4.10) 

where ap = — Xp/ \/l4. Over the range of Xp and parameter values considered, 
this curve differs httle from its small-Xp approximation, the hne Xs/Xv ~ 
-(k^^ - k'^^)/{k'^^ + + k'j^^). For any proton- neutron symmetric energy 
surface (k^^ = k^j,), the boundary curve (4.10) reduces to the hne xs = 0. 
The general problem of determining whether or not the prolate and oblate 
regions of the phase diagram share a common boundary involves solving for an 
intersection of the curves (4.7) and (4.10). For the proton-neutron symmetric 
case (k^^ = i^'vv)^ ^ prolate-oblate boundary arises for ^Ipp/ 1^'-^^, < (-\/l3— 2)/18 fa 
0.08919. 

The dependence of the boundary curve upon the relative strengths of like- 
nucleon and proton-neutron quadrupole coupling coefficients is shown in 
Fig. 10. As the n'^p are reduced in strength relative to k^^, the region of 
triaxial equilibrium contracts. Its separation from the origin and the onset of 
a direct prolate-oblate transition for k'^^/ k'^^^ < (vTS — 2)/18 is illustrated in 
Fig. 10(b). The parameter values chosen for Fig. 9(a) and for Fig. 10(a) are 
such that k'^j^/k'^^ is the same {n'^j^/n'^^, = 1/4) in both cases. Thus, though 
the curves in these two figures are quite different overall, by Eqn. (4.9) they 
share the same endpoint (x;/ ~ 0.6177) on the SU^i,(3)-SU*y(3) line. Numeri- 
cal results for the equilibrium values of the energy and coordinates along the 
SU7ry(3)-SU*^(3) line for different relative strengths of like-nucleon and proton- 
neutron quadrupole interactions are shown in Fig. 11. Although the transition 
to triaxial equilibrium is progressively delayed as k^^,/ k'^^, decreases, the triax- 
ial configuration of Fig. 1, with orthogonal symmetry axes {'Jt^ — O, ^y^ — Tr/S), 
is still ultimately achieved in the SU*^(3) limit. 

A qualitative understanding of the mechanism underlying the parameter de- 
pendences observed in Figs. 9 and 10 is easily obtained. Among the quadrupole 
terms in the energy surface, the (Qtt'Qtt) term has a dominant influence on the 
proton fluid equilibrium deformation, the {Q^, ■ Q^) terms drives the neutron 
fluid deformation, and the (Qw ■ Qu) term couples the two deformations. The 
SU*j,(3) conflguration with orthogonal symmetry axes (Fig. 1) arises since the 
{Qn ■ Qtt) term stabilizes the proton fluid about a prolate deformation and the 
{Qiy ■ Qv) term stabilizes the neutron fluid about an oblate deformation, while 
the (Qtt ■ Qv) term is responsible for the relative orientation of the symmetry 
axes. All three terms are essential to producing the SU*j,(3) triaxial equilib- 
rium configuration. If the strengths of the like-nucleon quadrupole interactions 
are both reduced relative to that of the proton-neutron interaction (Figs. 10 
and 11), the ability of the like-nucleon terms to stabilize the proton and neu- 
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Fig. 10. The curve of second-order phase transition between axial and triaxial de- 
formation in the plane ^' = 1, shown for different relative strengths of like-nuclcon 
and proton-neutron quadrupole interactions: (a) Htvk/k,-kv = K'wIK'-kv = 1/4 and 
(b) K-j^Tt / I'i-Ku = i^vv/ i^-Kv = 1/12. Both curves are for A^^r = The boundary curve 
for Ht^^/kt^i, = Khh/kt^i, = 1/2 (Fig. 3) is shown in each panel (dashed curve) for 
comparison. 

tron fluids about distinct deformations is reduced, against the tendency of the 
proton-neutron term to favor equal deformations. The onset of triaxiality is 
thus delayed. For Nj^/N^ greater than unity (Fig. 9), the relative weight of the 
{Qn ■ Qn) contribution to S increases, while that of the {Qu ■ Qu) contribution 
decreases. Thus, a larger larger positive Xu (oblate tendency for the neutron 
fluid) is needed for the neutrons to undergo a transition to an oblate config- 
uration against the restraining influence of the protons, while only a small 
positive Xn (oblate tendency for the proton fluid) is necessary for the protons 
to undergo a transition to an oblate configuration. 

Limited analytic results can also be obtained for the transition between spher- 
ical and deformed equilibrium. Consider the stability of the minimum at 
/^TT = /^i/ = 0. The second derivative of £ along the ray = and (3i, = 
is 

+ 2ul[el{l - eO - Kui'] - Su^u^KX cos(7. - 7.)- (4.11) 

As in the special case considered in Sec. 3.3, this quantity is independent of 
Xs s-nd Xv and depends upon the 7^ coordinates only through their difference 
77r — lu- As ^' is increased from 0, vanishing of this second derivative for some 
value of Ut,/ui, and ''i-K—lv indicates the onset of instability. Since the coefficient 
of cos(77r — 7;/) is negative, instability must first occur for 77r — 7i/ = 0. With 
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Fig. 11. Evolution of the equilibrium properties between the SU7rjy(3) and SU*j^(3) 
dynamical symmetries, for different relative strengths of like-nucleon and pro- 
ton-neutron quadrupole interactions: (left) k-j^t^/k-^i, = Ki^u/i^wv = 1/2 (for reference), 
(middle) k^^/k^j, = Kuu/nnu = 1/4 and (right) Kt^t^/k^^ = Ki,„/kj,i, = 1/12. The 
curves are for N^^ = iVjy. Proton fluid variables are represented by solid curves, while 
neutron fluid variables are represented by dashed curves. Points of phase transition 
are marked with arrows. 

Ttt — 7i/ set to zero, d'^£/dp'^\fs=o vanishes at 

Instability against deformations with equal and [3y therefore occurs at ^' = 
1/5 for any values of the energy surface parameters. However, the minimum 
at /9 = can in general become unstable against deformations with unequal 
/?7r and Pi, at a smaller value of The direction in (/S^r^/^i/) space in which 
instability first sets in is given by a quadratic equation for Ut^/u,^, 



^ 111/ ' 



(4.13) 



The ^' value at which instability occurs follows from this Uj^/ui, value by 
Eqn. (4.12). There are two important special cases in which instability first 
occurs at ^' — 1/5: (1) for the F-spin invariant Hamiltonian (3.2) with arbi- 
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Fig. 12. Evolution of the equilibrium properties along the U7rz/(5)-S07rz/(6) line in 

parameter space, illustrating instability against deformation with unequal /St^ and 
^y. Results are shown for the Hamiltonian (4.3) with e^^ = Sy and /Ctttt = i^vu = i^-ku/^ 
for (left) the symmetric case N^^/Ny = 1 and (right) the asymmetric case N^^/Ny = 4. 
The variable (3-,^ is represented by the solid curve, while fiy is represented by the 
dashed curve. The point of phase transition is marked with an arrow. 

trary N^^/Ny, and (2) for energy surface parameters which are proton-neutron 
symmetric {i.e., e'^ = and k'^^ = n'yy), as are obtained whenever a proton- 
neutron symmetric Hamiltonian is considered with NT^ = Ny. 

The evolution of equilibrium properties across the U^y(5)-S07r;/(6) second- 
order transition point for asymmetric in (3-,^ and [3y is shown in Fig. 12. 
From Eqns. (4.12) and (4.13), the phase transition for the parameter values 
used in the figure (e^r = ^vi i^tttt = i^uu = '^■nv/'^i and N^^/Ny = 4) occurs at 
^' ^ 0.1959, with P^/Py 0.7215. The evolution of the ratio P^/Py past this 
point is shown in Fig. 12(c). Observe that the proton and neutron equilibrium 
properties in this example are unequal even within the plane xv = 0, unlike 
the case of Sec. 3.1, illustrating that the the IBM-2 equilibrium problem in 
the general proton-neutron asymmetric case does not reduce to the one-fluid 
IBM problem. 

The minimum at zero deformation always becomes unstable at the ^' value 
determined by Eqns. (4.12) and (4.13), independent of and xv- However, 
this gives rise to a second order phase transition only if a first-order transition 
to a distinct minimum has not already occured at smaller ^' for those values of 
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Fig. 13. Phase diagram of the IBM-2 for the Hamiltonian (4.3), with 
i^tttt/i^-kv = i^vu/i^-Kv = 1/4 and Nt^/Nh = 4, as obtained by numerical minimiza- 
tion of £. The surfaces of first-order and second-order transition between regions 
of undeformed, axially symmetric deformed, and triaxially deformed equihbria are 
shown. The axes are defined as in Fig. 2. 

Xs and xvi as discussed in Sec. 3.3. For the F-spin invariant Hamiltonian (3.2) 
with arbitrary N^/N,^, an extension of relation (3.12) is readily found, namely 
that f(^', X5, xv; /3, 0, /3, 0) is invariant along any line of constant N^Xn + 
NyXv By arguments analogous to those of Sec. 3.3, a line of second order 
phase transition occurs for = 1/5 and X-kIXv = —N^/N^^, embedded in 
a surface of first order phase transition. For a Hamiltonian with arbitrary 
coupling constants, however, no such simple result is obtained. The full phase 
diagram of the IBM-2 Hamiltonian (4.3) for a set of parameters involving 
proton-neutron symmetric couplings (kj^tt/^^tti/ = i^w/ '^■kv = 1/4) but unequal 
boson numbers [N^/N^, — 4), obtained by numerical mimimization, is shown 
in Fig. 13. 

4-3 Majorana operator 

The Majorana operator M is the analogue in the IBM-2 Hamiltonian to the 

proton-neutron symmetry energy of the liquid drop model. This operator 
arises through a combination of direct shell model effects and renormaliza- 
tion effects [46,47]. As is evident from Eqn. (2.13), this operator energetically 
discourages configurations with ^ or ^T^^^y. The approximate strength 
of the Majorana contribution can be deduced from the energies of mixed sym- 
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metry excitations [48], including the scissors mode excitation in axial rotor 
nuclei {e.g., Ref. [49]). Comparison of the /^-vibrational or 7- vibrational en- 
ergy scale (~1 MeV) with the scissors mode energy scale (~2.5 MeV) indicates 
a; 5 to be a generally reasonable estimate. 

For axially symmetric configurations, the contribution of the Majorana oper- 
ator to the energy surface is, from Eqn. (2.13), (M) = Nt^Ni,{1 + f3'^)^'^{l + 
/5^)~^(/57r — PuY- If the equilibrium configuration without the Majorana op- 
erator already has /3,r — Pv-i as in the Xs — ^ plane for the proton- neutron 
symmetric case of Fig. 2, then introduction of the Majorana term has no fur- 
ther effect. Otherwise, the effect of the Majorana contribution is, naturally, 
to bring the equilibrium values /^^r and j3y closer to each other. In the plane 
^' = 1, this invalidates the simple results (3.6) and (3.7) giving the equilib- 
rium /Stt purely as a function of x-k and the equilibrium as a function of 
Xv The derivation of the simple equations (3.9) and (4.7) for the curve on 
which the axial minimum becomes unstable with respect to the 7p depended 
upon the use of these results to eliminate Xp iii favor of f3p at the minimum, a 
simplification which is no longer possible with a Majorana contribution. 

Investigation of the equilibrium properties in the presence of a Majorana op- 
erator must therefore rely upon numerical minimization. The evolution of the 
equilibrium properties between the SU^i,(3) and SU*j,(3) points in parame- 
ter space, for different strengths of the Majorana term in the Hamiltonian, 
is shown in Fig. 14. The Majorana operator is seen to have two main effects 
on the transition to triaxial equilibrium. The triaxial configuration of Fig. 1 
is highly proton-neutron asymmetric, and thus penalized energetically by the 
Majorana operator, as seen from Eqn. (2.13). Thus, first, the transition to 
towards such triaxial structure is delayed by the Majorana operator. Second, 
the Majorana operator has the effect of bringing the proton and neutron coor- 
dinate values closer together throughout the transition. The SU7r,y(3)-SU*j^(3) 
line is an extreme case. Without the Majorana operator [Fig. 5 or 14(a)], the 
evolution of the proton coordinates is completely decoupled from that of the 
neutron coordinates, and their equilibrium values are constant along the line. 
A Majorana strength of \/ k.^,^ ~ 10 is sufficent, however, to cause the two 
fiuids' (3p and 7^ coordinates to remain nearly equal throughout the evolution. 
The triaxial configurations produced are thus essentially one-fiuid, like those 
described by the Davydov model [50] or by the one-fiuid IBM with three-body 
or four-body operators in the Hamiltonian [26,28,51-53]. As the Majorana 
strength increases, the global minimum at triaxial deformation becomes very 
shallow, and the difference in energy between the ~ 7,^ ~ 7r/6 triaxial 
minimum and axial deformations (7^ = 7;/ = 0) approaches zero, as shown 
in Fig. 15. This indicates that the equilibrium structure for large Majorana 
strengths is essentially one-fluid ^-unstable, or S0(6)-hke, rather than rigidly 
triaxial. 
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Fig. 14. Evolution of the equiUbrium properties between the SU7r;/(3) and SU*^,(3) 
dynamical symmetries, for different Majorana operator strengths: (left) \/ k-j^v = 
(for reference), (middle) A/k^j^ = 1, and (right) = 10. The curves arc for 

the quadrupole coefficient values K^Tr/^Tn/ = Uvu/ i^-kv = 1/2 with N^^ = Ny. Proton 
fluid variables are represented by solid curves, while neutron fluid variables are 
represented by dashed curves. Points of phase transition are marked with arrows. 
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Fig. 15. Global minimum energy (solid curve) and the lowest energy for ax- 
ial (77r = 7i/ = 0) deformations (dashed curve), shown as a function of A/ Kf^jy , 
at the SU*^(3) point in parameter space. The difference is a measure of the 
7-stiffness of the triaxial minimum. The curves are for quadrupole coefficient values 

/^TTTr/^TTI/ = I^Vv/l^TTV = 1/2 Wlth iV^ = ■ 
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The decomposition of the Majorana operator into multipole components, eas- 
ily obtained from Ref. [21, (A2)], includes a nonzero hexadecapole contri- 
bution. The analysis of Ref. [21] therefore allows the existence of "oblique" 
equilibrium configurations, in which the proton and neutron intrinsic frames 
are not aligned. For axially symmetric aligned configurations {-di and 7^ all 
vanishing), instability with respect to the Euler angles is simple to investigate. 
The Hessian matrix of S with respect to the coordinates 7^, jjy, and 

■i?2 decomposes as a direct sum 7i = Tip © H-^ © Ti-^^^ with Tifj involving only 
derivatives with respect to the (3p, etc.. The dependence upon the azimuthal 
Euler angles i?! and ■i?3 vanishes by axial symmetry. Thus, instability would 
be indicated simply by a negative second derivative of S with respect to the 
opening angle between the proton and neutron frames. For the parame- 
ter ranges considered {k'^^ < 0, A' > 0, |xp| < \/7/2) and coordinate values 
encountered (/3p < \/2) in the present study, this quantity is always positive. 
The possibility of a first order transition to a distinct minimum with nonzero 
Euler angles must also be considered. Numerical searches provide no evidence 
for a such a transition, but such searches are necessarily not exhaustive due 
to the large number of parameters involved. 



5 Spectroscopic properties 

5. 1 Basic properties 

To provide a connection between the IBM-2 phase structure considered so 
far and observable quantities, let us now consider the evolution of the pre- 
dicted spectroscopic properties — eigenvalues and electromagnetic transition 
strengths — across the phase transitions. The present discussion emphasises 
the gross spectroscopic features which emerge in the transitions between the 
different regions of the phase diagram of Fig. 4. Numerical studies of several 
obscrvables in the vicinity of the SU*^(3) dynamical symmetry may also be 
found in Refs. [18,19]. 

Electric quadrupole and magnetic dipole transition matrix elements are cal- 
culated using transition operators [3] 



= eM + eA 



V 




and 




(5.2) 
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where Lp = \^{dj, x dp)'^^\ The transition strengths are B{aX; Ji — > J/) = 
(2 J/ + 1) |(/||T('^^)||i)|V(2Ji + 1). Schematic values = = 0.1 eb, = 
/iAT, and (7jy = are used for the effective charges (see Ref. [9] for further 
discussion). The Xp"^"^^ parameters are taken equal to their counterparts Xp 
in the Hamiltonian, following the consistent quadrupole formalism [54]. The 
Hamiltonian used for the present discussion is the simple form (3.2) with the 
addition of the Majorana operator. Diagonalization is carried out using the 
computer code npbos [55], for boson numbers — Ni, — 5. 

The energy spectrum for the SU*j,(3) symmetry is known analytically [10, 17]. 
The operator (Q^^ +0^") ■ (Q^^+Qi^) with = -V7/2 and Xu = +V7/2 can 
be reexpressed in terms of the quadratic Casimir operators [3] of the subalgebra 
chain (2.5). as (3/4)^2 [SU;, (3)] - (3/16)C2[SO^,(3)]. For H = -Q ■ Q, the 
eigenvalues are thus given by 

E{X, At, L) = -^(A' + f/ + Xfx + 3X + 3//) + |l(L + 1), (5.3) 
2 o 

where (A,/i) are the SU*j,(3) Elliott quantum numbers [56] and L is the 
angular momentum quantum number. This is the SU(3) energy formula, 
and level energies within an SU*j,(3) representation follow L{L + 1) rota- 
tional spacings. However, the allowed SU(3) representations (see Refs. [10, 17]) 
are different from those for 811^,^(3). The ground state representation is the 
{2Nj^, 2Niy) representation, rather than the usual 811^,^(3) ground state repre- 
sentation {2N^+2N^, 0). According to the SU(3) D SO (3) branching rules [56], 
the ground state representation thus contains multiple degenerate rotational 
bands, with K quantum numbers 0, 2, 4, . . ., as shown in Fig. 16. [It has be- 
come customary to label the bands of an SU(3) representation by the Elliott K 
quantum number [56], in part by analogy with the rigid rotor spin-projection 
quantum number K [35], which yields a band of the same band head spin. 
However, the Elliott SU(3) basis is not orthogonal, and the orthonormal states 
obtained from diagonalization of a Hamiltonian are thus not in general Elliott 
states. Throughout this article, the label K is used simply to indicate the spin 
content of a band.] Other SU*^(3) representations appearing at low energy 
include the (2iV^ - 4, 2N^ + 2), {2N^ + 2, 2N^ - 4), and (2iV^ - 1, 2iV^ - 1) rep- 
resentations (Fig. 16). In the classical limit, the first two of these correspond 
to coupled P and 7 vibrations of the fluids [20]. The last corresponds to an 
"orthogonal" scissors mode, in which the proton and neutron symmetry axes 
oscillate about their equilibrium perpendicular relative orientation [17,20]. 

The B{E2) strengths for transitions between levels of the SU*j,(3) ground 
state representation, calculated numerically, are shown in Fig. 17(a). These 
strengths follow selection rules dictated by the presence of a discrete, parity- 
like symmetry. Consider the operation consisting of negation of the dp bosons 
followed by interchange of all proton and neutron bosons, together yielding 
<-> Si, and <-> —d,^. This was denoted the T parity operation by Ot- 
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Fig. 16. Level scheme for the SU*j^(3) dynamical symmetry, following the energy re- 
lation (5.3). The lowest-energy 8X1*^(3) representations, labeled by the (A,//) quan- 
tum numbers, are shown for Nt^ = Ni, = 5. 

suka [57]. For 

Xtt — Xvi the IBM-2 Hamiltonian (2.1) with e-^ — Sy and 
is invariant under the T parity operation. Thus T parity is a good symmetry 
throughout the central vertical plane of Fig. 2, including the SU*^(3) dynam- 
ical symmetry. The operator T*^^^^ carries negative T parity. Therefore, E2 
transitions occur only between states of opposite T parity, and all electric 
quadrupole moments vanish. The operator t'^^^) decomposes into a part of 
positive T parity oc (Ltt-I-Lj,), which generates magnetic dipole moments, and 
a part of negative T parity oc (L^r — I/^), which induces transitions between 
different states [57]. Thus Ml transitions follow the same T parity selection 
rule as E2 transitions, but magnetic moments are allowed. In applying the 
selection rules, it should be noted that the E2 selection rule only holds in its 
exact form for e.,^ — e^. Also, the T parity operation involves interchange of 
all proton and neutron bosons, so it is only well defined for — N^, but the 
selection rules persist approximately even for N^^^N^. 

The classical interpretation of the T parity operation is obtained by observing 
that, in the definition of the coherent state (2.7), negation of the dp boson is 
equivalent to negation of the a^^j^ deformation coordinates. In the geometric 
model, this is the 7 parity operation of Bes [58], which exchanges prolate 
and oblate liquid drop deformations. Thus, the T parity operation exchanges 
prolate and oblate deformations for each fluid and then interchanges proton 
and neutron fluids. The SU*^(3) triaxial conflguration of Fig. 1 is seen to be 
invariant under this combined transformation. 

We now note the general characteristics of the B{E2) strengths for the SU*j,(3) 
symmetry shown in Fig. 17(a). While the energy levels for the 811*^,(3) sym- 
metry fall naturally into rotational quasi-bands, the transition strengths are 
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Fig. 17. Level schemes and selected B{E2) strengths for (a) the SU*j,(3) dynamical 
symmetry {xn = -a/7/2 and Xu = +a/7/2), (b) X/n.^^ = 1, and (c) X/k^^u = 10. 
The level scheme for the Davydov model with 7 = 7r/6 is shown for comparison (d). 
Levels of positive and negative T parity are indicated by solid and dashed lines, 
respectively. All B{E2) strengths are given for the spin-descending transition di- 
rection, to allow direct comparison between panels, regardless of the ordering of 
level energies, and are normalized to B{E2;2'l 0^). Energies are normalized 
separately in each panel. A minimal Majorana term {X/ K-j^y = 0.125) is used in the 
calculation of part (a) to lift the degeneracy of the SU*^(3) multiplets. Calculations 
are for AT- = Ny = 5. 
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far from those for an SU,ri/(3) rotor. Many of the interband transitions have 
strenths of the same order as in-band transitions, while some of the in-band 
transitions vanish due to the T parity selection rule. (Levels of positive and 
negative T parity are indicated by sohd and dashed hnes, respectively, in 
Fig. 17.) Transitions between non-adjacent bands, not shown in the figure, are 
weaker by an order of magnitude or more. 

Magnetic dipole transitions arise in collective models from separation be- 
tween the proton and neutron fluid distributions {e.g, Ref. [59]). The static 
asymmetry between these distributions for the SU*^(3) dynamical symme- 
try gives rise to extremely large Ml admixtures whenever spin-allowed, with 
B{M1)/ B{E2) ~ l/x^/(e^b^) for most of the transitions in the ground state 
representation. These Ml strengths are comparable to the Ml decay strength 
of the scissors excitation. 

The SU*^(3) dynamical symmetry has previously been loosely associated with 
the 7 = 7r/6 triaxial rotor of the Davydov model [50], on the basis of two 
similarities: the angular momentum content of the ground state representation 
and the quadrupole moment components found in the classical limit [10]. The 
level energies in the two cases differ considerably, as the Davydov model bands 
exhibit a characteristic (2+3+) (4^5"'") . . . clustering or "staggering" of energies 
[Fig. 17(d)]. However, a detailed comparison of the E2 transition strengths of 
the two models [Fig. 17(a,d)] reveals an extraordinary similarity, with many 
E2 strengths differing by only a few per cent. The full relationship between the 
models has not been established (an approximate SU(3) algebra underlying 
the dynamics of the triaxial rotor has been discussed in Ref. [60] ) . 

If the Majorana operator is introduced into the SU*^(3) Hamiltonian, the en- 
ergy spectrum is radically altered, as depicted in Fig. 17(b,c). The F-spin [8], 
formally analogous to isospin, of a state provides a measure of its proton- 
neutron asymmetry. For sufficiently pure F-spin F = F^ax = ^{^tt + N^), 
the IBM-2 effectively reduces to the one-fluid IBMJ61]. The SU;^(3) eigen- 
states are highly mixed in their F-spin content ((F^)^/^ O.SFmax), reflect- 
ing the proton-neutron asymmetry of the classical equilibrium configuration. 
The Majorana operator, by energetically penalizing states of F < Fmax, pu- 
rifies the F-spin contents of the low-energy states. E.g., the SU*^(3) ground 
state is left with an F < F^ax content of only about ~ 1% for \/ n-j^i, = 10 
[Fig. 17(c)]. A small Majorana contribution in the Hamiltonian, X/kj^-i^ < 1, 
results in a (2"''3"'")(4"'"5+) . . . staggering of level energies [Fig. 17(b)] like that 
of the Davydov model [Fig. 17(d)], as noted in Ref. [18]. A larger Majo- 
rana contribution [Fig. 17(c)] leads to the reverse 2"'"(3"'"4"'")(5"'"6"'") . . . stag- 
gering, characteristic of 7-soft structure. The F2 strengths evolve toward 
S0(6)-hke values [3] with increasing Majorana strength, and the Ml ad- 
mixtures decrease rapidly from their large SU*^(3) values, with a typical 
scale B{M1)/B{E2) ~ 10-Vir/(eV) for Fig. 17(b) and B{M1)/B{E2) ~ 
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Fig. 18. Excitation energies of the lowest 2+ (solid), 3+ (dashed), and 4+ (dot- 
ted) levels along the SU7ri/(3)-SU*j,(3) transition, (a) with no Majorana operator 
and (b) for X/kt^u = 5. The change in slope for the excited band level energies at 
~ —0.4 in part (a) arises from the crossing of the scissors and 7 bands. Calcula- 
tions are for x-k = — and = = 5, with minimal perturbations ^' = 0.95 
and for part (a) A/ktt!/ = 0.125 to remove numerical degeneracies. 

10~^/i|f/ (e^b^) for Fig. 17(c). The evolution observed of the apparent structure 
with increasing Majorana contribution — from SU*j,(3) triaxialiaty, through 
one-fluid rigid triaxiality, to one-fluid 7-softness — is as expected from the 
classical limit analysis of Sec. 4.3. 

5.2 The ^\5t,p{?>)-^\51^{?>) transition 

The evolution of the IBM-2 predictions along the SU7ri,(3)-SU*^(3) transition 
is shown in Figs. 18-20. Since the extremely low lying K — 2 band is char- 
acteristic of the SU*^(3) level scheme, the evolution of the relevant energies 
is plotted in Fig. 18. For simplicity, we first consider the transition without a 
Majorana interaction. In this case, rotational L{L+ 1) energy spacings within 
bands are almost exactly preserved throughout the transition. 

In the SU7rj/(3) limit [Fig. 19(a)], the lowest excitation above the ground state 
band {F — Fj^ax) is the scissors mode {F — Fmax — 1), followed at higher energy 
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Fig. 19. Level schemes for the SU,r!/(3)-SU*j^(3) transition with no Majorana oper- 
ator: (a) the SU7r;y(3) hmit at = — '\/7/2, (b) the second order transition point at 
Xu ~ 0.4035, (c) the approximate ground state crossing point at Xf = 0-8, and (d) the 
SU*j,(3) Hmit at Xf = +'\/7/2. The E2 branching patterns are shown for selected 
band head states, and the ground state band 4^ 2"*" and 2"*" O"*" transitions 
are included for scale [shaded arrows, widths proportional to B{E2) strength]. The 
Ml transitions deexciting the scissors band head states are shown as well [dashed 
arrows, widths proportional to B{M1) strength]. The g (ground), /?, 7, and 6 (scis- 
sors) band designations are schematic (see text). All panels share a common energy 
scale. Calculations are for Xn = — \/7/2 and = Ny = 5, with minimal pertur- 
bations ^' = 0.95 for part (a) and X/k-j^v = 0.125 for part (d) to remove numerical 
degeneracies. 
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Fig. 20. Level schemes for the SU7ri/(3)-SU*j^(3) transition with Majorana contribu- 
tion A/ktti/ = 5: (a) the SU7ri/(3) hmit at Xv = —V^/"^, (b) the second order transition 
point at Xiy ~ 1-022, and (c) the SU*j^(3) hmit at Xf = +\/7/2. For the transition 
point, B{E2) and i?(Ml) strengths among selected low-lying levels are given if 
part (d). Parts (a)-(c) share a common energy scale, while energies in part (d) 
are normalized to the 2^ level energy. Dashed arrows in parts (a)-(c) indicate Ml 
transitions from the scissors excitations. B{E2) and B{M1) strengths in part (d) 
are shown spin-descending, with B{E2) strengths normalized to B{E2;2^ 0^) 
[0.283 e^b^]. The g (ground), /5, 7, and 6 (scissors) band designations are schematic. 
Calculations are for x-k = —V7 /2 and N.,^ = Ni, = 5, with a minimal fid perturbation 
^' = 0.95 for part (a) to remove numerical degeneracies. 

by a cluster of degenerate K — {) and 2 excitations. These arise from proton- 
neutron symmetric beta and gamma oscillations {F = Fmax), asymmetric beta 
and gamma oscillations {F = F^^^x — 1), and two-phonon scissors oscillations 
(F — Fmax — 2). All E2 transitions from the excited bands to the ground 
state band vanish exactly for the SU7ri/(3) dynamical symmetry, but for small 
breakings of the symmetry they approximately follow the rotational Alaga 



36 



rules. Ml transitions vanish as well, except those involving excitations of the 
scissors mode. 

The level scheme for the SU7ri/(3)-SU*j,(3) second order transition point {xu ~ 
0.4035) is shown in Fig. 19(b). As Xi^ is increased from — \/7/2, one of the 7 
bands (no longer of pure F spin) rapidly descends to lower excitation energy. 
At the transition point, this band is connected to the ground state band by 
E2 transitions ~ 1/10 as strong as in-band transitions. The Ml admixtures 
are extremely large, with B{M1)/B{E2) ~ l/x^/(e2b^), both for the 7 to 
ground and the in-band transitions. The 7 band is followed in its descent by a 
series of K = 4, 6, . . . bands, originating as two-phonon, threc-phonon, etc., 7 
excitations. However, the K = two-phonon 7 excitation docs not descend as 
rapidly and thus has a positive anharmonicity, ~ 2.8 at the transition point. 
Positive anharmonicity of this band has also recently been discussed as a 
signature of the axial-triaxial transition both in the one-fluid IBM [62] and in 
the Y(5) geometric model [63]. [The irregular level spacing seen in this band 
in Fig. 19(b) is not a fundamental signature but rather arises from mixing 
with the nearly degenerate scissors band for these parameter values.] As the 7 
vibrational energy scale decreases, the K — 1 and 3 7 exitations of the scissors 
band similarly approach the scissors band in energy. 

Although the coherent state energy surface first becomes soft to triaxial defor- 
mations at the second order transition point, it is only for X;/ > 0.65 that highly 
triaxial configurations like that of Fig. 1 become lower in energy than axial 
configurations. A level scheme obtained in this regime, at Xu = 0.8, is shown 
in Fig. 19(c). The ground state band and first K = 2 and 4 bands are close 
in energy, resembling the SU*j,(3) ground state representation. The relative 
E2 transition strengths are intermediate between those for the SU^j^(3) and 
SU*j,(3) symmetries. The first excited K — band and following K — 2 band, 
in contrast, have relative energies more appropriate to a ground state and its 
associated 7 excitation, and B{E2) values closely matching the SU7rj/(3) Alaga 
rules. The energy spacing scales within these bands is less (by ~ 10%) than 
that of the ground state family of bands, suggesting a different moment of 
inertia or different triaxial nature to the rotation. The excitation energy of 
the O2 state reaches its lowest value along the SU,ri/(3)-SU*y(3) transition at 
approximately this Xu value. The observables are thus highly suggestive of an 
avoided crossing between an SU*^(3)-like configuration and an SU7ri/(3)-like 
configuration, although the actual eigenstate structure is presumably more 
complicated. An apparent avoided crossing also occurs between 1+ scissors 
excitations based upon these configurations. 

We now briefly reconsider the SU^jy(3)-SU*^(3) transition for a realistic Ma- 
jorana operator strength (A/k^j^ = 5) [Figs. 18(b) and 20]. At the SU^i.(3) 
limit, the eigcnstatcs are unchanged, but the levels with F < F^^x are raised 
in energy [Fig. 20(a)]. The lowest-lying excitations are thus proton-neutron 
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symmetric (3 and 7 excitations, followed at higher energies by the symmetric 
P and 7 multiphonon excitations. The second-order transition between axial 
and triaxial equilibrium occurs at much larger Xu than without the Majo- 
rana operator, Xu ~ 1.022 (Fig. 14). The level scheme for the transition point 
[Fig. 20(b)] shows large differences between ground state and excited famihes 
of levels. The lowest K = 0, 2, and 4 quasi-bands have level spacings inter- 
mediate between those for the SU7rjy(3) and SO.,ru{Q) symmetries. The next 
K = quasi-band is far less rotational in spacing and has a 2+-0"'" energy 
spacing about twice as large as for the ground state. The levels above it are 
assembled into an approximate S07ri/(6) multiplet structure. All Ml strengths 
are attenuated relative to the cases with no Majorana contribution, with a 
typical admixture scale B{M1)/ B{E2) ~ 0.05;Li^/(e^b^). Qualitatively, little 
further changes as Xu increases to +\/7/2 [Fig. 20(c)], except that the ground 
state family of levels takes on more distinctly S0^,y(6)-like energy spacings 
and B{E2) strengths. 



5.3 The U7ri/(5)-SU*y(3) transition 



The basic spectroscopic features along the U7r;y(5)-SU*j,(3) transition differ 
fundamentally depending upon whether or not a realistic Majorana contribu- 
tion is included in the Hamiltonian. The level scheme for the U7riy(5) symmetry 
of the IBM-2 [3, 9] exhibits extreme degeneracies if no Majorana operator is 
present: there are two degenerate one-phonon 2"*" levels (one fully symmetric 
and one of mixed symmetry) and eleven degenerate levels of various symmetry 
in the two-phonon multiplet. By the second order transition point ((^' = 1/5), 
these highly degenerate levels evolve into a cluster of low-lying excitations of 
widely varying F-spin content, shown in the inset to Fig. 21(a). Among these 
levels, incipient forms of the 811*^(3) ground-state and scissors quasi-band 
structures are apparent. The electromagnetic transition properties are dom- 
inated by the ubiquity of strong Ml transitions, arising from the mixing of 
U7r!/(5) states of different F spin. For instance, both the 2 J — > 2^ and 23 — > 2^ 
transitions have B{M1) strengths > 0.1 typical of mixed symmetry state 
decays. 

If instead a realistic Majorana operator strength \/ n-r^y ~ 5 is considered, the 
low-energy excitation spectrum simplifies considerably, as shown in Fig. 21(a). 
The low- lying levels have nearly pure F Fmax, with a few per cent admix- 
tures of other F values, and are grouped into approximate 5*0(5) multiplets. 
Mixed symmetry levels oi F ^ i^max — 1 ai'e present at higher energy (1"*", 
2"*", and 3"*" levels at top of figure). The T parity is a good quantum number 
throughout the U7ry(5)-SU*j,(3) transition, so T parity selection rules apply to 
the electromagnetic transitions. The level scheme and E2 transition strengths 
are nearly identical to those obtained for the second order transition point 
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Fig. 21. Level schemes for the second order transition points along (a) the 

U7r!/(5)-SU*^(3) transition and (b) the U7r;y(5)-S07r,y(6) transition, calculated for 
a Majorana operator strength of A/k^j^ = 5. The inset to part (a) shows the corre- 
sponding level scheme with no Majorana operator. Energies are indicated normal- 
ized to the 2'^ level energy. B{E2) strengths are shown spin-descending, normalized 
to B{E2;2f Of) separately in each panel [0.175 e^b^ for part (a), 0.168 e^b^ 
for part (b)], and selected B{M1) strengths are indicated as well. Dashed double 
arrows indicate selected strong Ml transitions. Calculations are for A^^ = A^^ = 5, 
with Xtt = -V7/2, Xu = +-\/7/2, and ^' = 0.2 for part (a) and with = Q, = 0, 
and C' = 0.2 for part (b). 
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between U^,y(5) and S07ri/(6), shown for comparison in Fig. 21(b). There is 
thus also considerable resemblance to the predictions of the geometric E{5) 
model [64,65]. Two features of the U7ry(5)-SU*j,(3) transition point spectrum 
allow it to be distinguished from the U^,^(5)-SO^,^(6) case. First, a slight break- 
ing of the degeneracy of the SO (5) multiplets occurs, in the sense necessary 
to create the SU*j,(3) quasi-bands. Thus, the 22 level moves below the 4f 
level, the 3f and 4^ levels move below the 6^ level, etc. This effect becomes 
much more marked just past the transition point. Second, sizeable Ml transi- 
tion strengths, B{M1)/ B{E2) ^ 0.05 i^%/{e'^h^) to 0.1 i^%/{e^h^), are present 
among the low- lying levels [Fig. 21(a)], while these vanish when Xs — ^- Such 
Ml admixtures in general play an essential role in identifying deviations from 
XV ^0 [57,66,67]. 



6 Conclusion 

The IBM-2 phase diagram investigated here provides a framework for study- 
ing the transition between axial and triaxial structure in nuclei. Triaxial de- 
formation might arise from several possible sources: distinct deformations of 
the proton and neutron fluids as considered in the present work, higher-order 
interactions in an essentially one- fluid nucleus [28], or the presence of conflg- 
urations involving hexadecapole nucleon pairs [68]. The main spectroscopic 
features of proton- neutron triaxial structure are exemplified by the SU*j,(3) 
dynamical symmetry, shown in Fig. 17(a). Proton-neutron triaxiality is char- 
acterized by a low-lying K — 2 band, as in the other forms of nuclear triaxiality, 
but with level energies following a rotational L[L + 1) sequence. The B{E2) 
strength pattern is remarkably similar to that of the classic rigid triaxial rotor 
of the Davydov model, including the unusual feature that B{E2] 22 Of) = 
(this however only strictly holds for A^^ = N,y and e^, = e^). The microscopic 
conditions leading to two-fluid triaxial structure, which requires X'k and Xv 
of opposite sign, are found when the proton bosons are particle like (below 
mid-shell) and the neutron bosons are hole-like (above mid-shell) or vice versa. 
Nuclei in this category include the heavy Ru and Pd isotopes [67,69, 70], the 
light rare earth nuclei below the N — S2 shell closure, and the extremely light 
W, Os, and Pt nuclei. 

The SU*j,(3) dynamical symmetry probably does not occur in its pure form 
in any actual nuclei. In this article, we have seen how the features of 
SU*^(3) proton- neutron triaxiality are modifled by the Majorana interaction 
[Fig. 17(b,c)], which tends to alter the spectrum towards that for the S0^,^(6) 
dynamical symmetry. Of most importance from a phenomenological viewpoint, 
we have investigated the transition from axially symmetric deformed [SU7r;y(3)- 
like] to triaxially deformed [SU*y(3)-like] structure, either without [Fig. 19] or 
with [Fig. 20] a Majorana interaction. In particular. Fig. 20 represents a re- 
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alistic study of this transition, including the behavior of the lowest K = 1, or 
scissors, band. It is found that, despite the attenuation of proton-neutron triax- 
ial features by the Majorana interaction, two-fluid triaxiality can serve as the 
underlying mechanism driving one-fluid triaxiality, and recognizable features 
of its two-fluid nature can remain even in the presence of realistic Majorana 
interactions. The two-fluid phenomena are likely to take on more importance 
as very neutron rich nuclei become accessible to experimental study. When 
the valence protons and neutrons occupy well-separated orbitals, the proton- 
neutron interaction strengths are reduced, yielding a situation seen in Sec. 4 
to be much more conducive to true two-fluid triaxiality. 

The present analysis may also serve as a model for the study of other multi- 
fluid bosonic systems. Within nuclear physics, the Ucore(6) ® Uskin(6) descrip- 
tion of core-skin collective modes in neutron rich nuclei [71] is most directly 
analogous. In molecular physics, the C/i(4) f/2(4) vibron model with two vi- 
bronic species [4] may be studied similarly. The appropriate coherent state 
formalism is discussed in Ref. [72]. The phase transition in this model is 
relevant to coupled vibronic bending modes in "floppy" molecules such as 
acetylene [73]. An extreme case of a multi-fluid algebraic model is found in 
applications to polymers, in which an arbitarily large number of vibronic flu- 
ids are coupled [74]. Yet another case is that of atomic condensates, for which 
scissors modes of a single-constituent Bose-Einstein condensate relative to an 
anisotropic potential have been observed [75] . Experiments have been planned 
to produce condensates of two atomic species. The exotic features of the IBM- 
2 phase diagram are likely to be encountered for other multi-fluid systems as 
well. They highlight the need for a classification scheme beyond the simple 
Ehrenfest or one-parameter Landau models to address the phase structure of 
systems which simultaneously possess multiple control parameters and multi- 
ple order parameters. 
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A Matrix element of an cirbitrary m-body operator between multi- 
fluid coherent states 



For calculations in the coherent state formalism involving multi-fluid oper- 
ators, it is useful to deduce a general formula for the matrix element of an 
arbitrary m-body operator between arbitrary multi-fluid coherent states. De- 
flne S different coherent boson species Bl {s = 1, . . . , S) in terms of the basic 
boson operators bl {i — l,...,n) as orthonormal linear combinations 



Bl 



n 

1=1 



(A.l) 



where a 



S' ,1 



a 



s,i — Sg's- Then the multi-species coherent states 



|A^i---A^5> 



s 

n 

s=l 



{BirA\o) 



(A.2) 



arc normalized, are orthogonal to each other for different A^^ values, and have 
good total boson number = l]f=i Ng, where the total boson number operator 
is A" = J2i=i blh. Coherent states of the type (A.2) also arise in the study of 
intrinsic excitations in algebraic models. In this context, the different coherent 
boson operators Bl represent a ground state condensate and one or more 
orthogonal excitation modes (see, e.g., Refs. [72,76-79]). 

The matrix element of an arbitrary m-body operator (m > 1) between two 
arbitrary multi-species coherent states is [80] 



{K---N's\[m'^)[EK)\Ni---Ns) 

S r S I — — ; -I , m X 

E n<^^i-^iv.-..m^ArJ^ n«tl,.^«*.n , (A.3) 

t' t' , '-5=1 -' \=\ ^ 



r f 

where Vs = J2^iSu,s, = J2^i^t'.,s, and an underlined superscript indicates 
the falling factorial [m- = m(m — 1) • • • (m — n -|- 1)] [81]. The multiple sum 
in (A.3) nominally contains S'^"^ terms, but only those with Nl — iy'^^Ng — Vg 
and Ng — i/s > can yield nonzero contributions. Three stages are involved 
in evaluating the matrix element of a general operator: reexpression of the 
operator in terms of normal-ordered m-body terms, evaluation of the matrix 
elements of these by (A.3), and simphflcation of the result. These steps can 
all readily be carried out though computer-based symbolic manipulation. 
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B Properties of F{a; x) 



The IBM-2 energy surface considered in Sees. 3 and 4 involves terms of the 
form 

_ a;(l + ax) 
1 + x^ 

We summarize here the extremum properties of this expression. For any value 
of the parameter a, F{a;x) has two extrema with respect to at 



Fia;x)^--^^^. (B.l) 



x±{a) = a±VaFTl, (B.2) 

The global minimum of F{a; x) is located at a;_(a), and the global maximum 
is located at x^{a). Note that x_{a) < 0, x^{a) > 0, and the two extremum 
positions are related by X-{a)x+{a) = —1. The extremal values of F are simply 

F[a;x±{a)] = ^x±{a). (B.3) 

If a > then |F[a; x_(a)] | < |F[a; x+(a)]|, while for a < the inequality holds 
in the opposite sense. 
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